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SUMMARY 


A theoretical model developed to calculate pressure-temperature cross 
spectra, pressure spectra, temperature spectra and pressure cross spectra In a 
ducted combustion system Is presented- The model assumes the presence of a 
f luctuatlng-volumetrlc-heat-release-rate disk source and takes Into account 
the spatial distribution of the steady-state volumetric-heat flux. Using the 
model, pressure, velocity, and temperature perturbation relationships can be 
obtained. The theoretical results show that, at a given air mass flow rate, 
the calculated pressure-temperature cross spectra phase angle at the combustor 
exit depends on the model selected for the steady-state volumetric-heat flux 
In the combustor. Using measurements of the phase angle, an appropriate source 
region model was selected. The model calculations are compared with the data. 
The comparison shows good agreement and Indicates that with the use of this 
model the pressure-temperature cross spectra measurements provide useful In- 
formation on the physical mechanisms active at the combustion noise source. 
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NOMENCLATURE 

see Eq. (46) 
area, 

see Eq. (30) 
see Eq. (47) 
see Eq. (45) 

sound propagation speed, m/s 

gas specific heat at constant pressure, J/kg-K 

gas specific heat at constant volume, J/kg-K 

diameter ,m 

frequency, Hz 

see Eq. (25) 

see Eq. (24) 

<-D 1/2 
see Eq. (54) 
summation Index 
see Eq. (26) 
see Eq. (26) 
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see Eq. (27) 
distance between points 
summation Index 
see Eq. (28) 

Mach number 

gas molecular wt: 

number of regions 

see Eq. (45) 

acoustic pressure. Pa 

volumetrlc-heat-rel ease-rate coefficients 

volumetric-heat-release rate 

gas constant 

pressure reflection coefficient 

nozzle exit radius 

see Eq. (54) 

gas entropy, J/kg-K 

transfer matrix 

temperature, K 

mean flow gas velocity 

mass flow rate, kg/s 

cartesian coordinate 

see Eq. (29) 

acoustic Impedance 

specific heat ratio of gas 

unit Impulse function 

see Eq. (56) 

see Eq. (62) 

see Eq. (65) 

time, s 

see Eq. (51) 

see Eq. (55) 

gas density, kg/m 

see Eq. (65) 

angular frequency, 2*f 


V 

2 



complex conjugate 
transpose 

Instantaneous quantity 


reference quantity 
perturbed quantity 
combustor Inlet 
duct exit 
gas 

duct Inlet 
source 


INTRODUCTION 

Combustion noise occurs In many systems (refs. 1 and 2) and Is an Impor- 
tant noise source In turbine-powered aircraft and helicopters (refs. 3 and 4). 
The NASA Lewis Research Center turbine engine core noise research program In- 
volved tests In turbine engines, engine component facilities, and small scale 
rigs using In-house, contractor, and university facilities. As part of this 
research program, various measurement techniques were developed to study noise 
associated with combustion rigs and aircraft engines. The experimental tech- 
niques using unsteady pressure measurements to study core noise Included: co- 
herence measurements (refs. 5 and 6), conditional spectral analysis (ref. 6), 
modal analysis (refs. 7 and 8), three signal coherence analysis (refs. 9 to 
11), cepstral analysis (ref. 12) and pressure-temperature cross spectra 
(refs. 13 and 14). Note that the last technique also Includes temperarture 
measurements. 

In analyzing combustion noise measurements two common assumptions are 
made: the same basic source mechanism Is present for all tests, and the mag- 

nitude of the noise depends on a single parameter, related to the heat-release 
rate, that can be calculated from the steady state operating conditions and 
the geometry. Thus, the existence of a relationship between the combustion 
noise and the combustion process has generally been established by determining 
that a correlation between the combustion noise and a steady-state heat-release 
rate parameter exists (refs. 10, 12, and 15 to 17) rather than by any direct ' 
physical measurement of a fluctuating quantity. The parameter used to corre- 
late combustion noise can be selected using a theoretical model (ref. 17) or 
Judicious curve fitting (refs. 10 and 16). In either case, the effectiveness 
Is decided by examining the ability of the parameter to correlate noise meas- 
urements at many test conditions. 

Unfortunately, this approach does not yield any Information about possible 
physical mechanisms active at the combustion source. Note that this Is In 
strong distinction to the case of flow noise due to a solid body In contact 
with airflow containing velocity and pressure fluctuations. In this case, the 
various physical mechanisms are known and the assumption Is made that predict- 
ing flow noise requires the specification of air stream mean velocity, rms 
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turbulence Intensity, correlation lengths, and turbulence spectra (ref. 18). 
However, most of these quantities can not be determined from the mean operating 
conditions and test geometry but have to be measured at each test condition. 

\ ’ . J 

The Importance of a physical mechanism Is tested by Including It In a 
theoretical model which Is used to calculate results that can be compared with 
measurements. Analytical one-dimensional models for combustion noise genera- 
tion and propagation In constant and variable area ducts were developed both 
at NASA Lewis and with NASA Lewis support (refs. 19 to 28) to compare with 
experimental measurements. Even though these models yield acoustic pressure 
spectra, cross spectra and transfer functions that resemble measured data this 
agreement Is not sufficient to determine which source region model Is appro- 
priate since similar acoustic results are obtained for different source region 
models. 

Measurements were made of pressure-temperature cross spectra (refs. 13 to 
14) along with pressure auto spectra and cross spectra to select an appropriate 
source model and to provide Information on the physical mechanisms In the com- 
bustion noise source region. The simultaneous real-time measurement of tem- 
perature and pressure fluctuations In a duct Is not new or unique and neither 
Is the derivation of auto spectra from them (refs. 29 to 32). However, the 
use of the cross spectra between pressure and temperature to show that a phys- 
ical connection between the heat-release rate and the combustion noise does 
exist Is new and unique (refs. 13 and 14) as Is the use of this data to develop 
a model. 

In the next section, the model developed to analyze measurements of 
pressure-temperature cross spectra, pressure spectra, temperature spectra and 
pressure cross spectra made In a combustion duct Is presented. Then the ex- 
perimental Investigation Is discussed and the model results are compared with 
experimental measurements. 


ANALYTICAL MODEL 

A theoretical model developed to calculate pressure-temperature cross 
spectra, pressure spectra, temperature spectra and pressure cross spectra In a 
ducted combustion system Is presented. The model assumes the presence of a 
fluctuatlng-volumetrlc-heat-release-rate disk source and takes Into account 
the spatial distribution of the steady-state volumetric-heat flux. A solution 
Is obtained by solving an Inhomogeneous differential equation using three 
boundary equations. One boundary condition Is obtained by assuming the Inlet 
gas behaves adlabatlcally. A second Is based on using a pressure reflection 
coefficient at the combustor entrance. The third Is due to assuming a flanged 
pipe pressure boundary condition at the exit. Using the model, pressure, 
velocity, and temperature perturbation relationships can be obtained. The 
model equations used are similar to those discussed In references 20 and 21. 
The model developed to calculate pressure-temperature cross spectra Includes 
terms dependent on a space-dependent-volumetric-heat-source region In addition 
to terms dependent on area gradients. While, the experimental data were ob- 
tained In a constant area combustion rig, the equations presented are for the 
more general case. In references 20 and 21, these equations were used only to 
calculate pressure spectra. In this paper, they are used to generate tempera- 
ture spectra and pressure-temperature cross spectra. The equations are solved 
using a transfer matrix approach rather than by direct numerical Integration 
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as In references 20 and 21. This method gives perhaps more physical Insight 
at the expense of some additional analytical work. However, Including a dis- 
tributed volumetric heat-release rate perturbation source region Is also more 
difficult. Consequently, the calculations are for a single disk source 
volumetric-heat-flux-rate perturbation. Models have been developed previously 
for sound propagation In a duct with an axial temperature gradient (refs. 33 
to 39), however, they do not Include a heat source. 


Governing Equations 


The one dimensional continuity, momentum, and energy equations for the 
bulk gas are: 
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where $, p, 0, t, and q are respectively the Instantaneous pressure, 
density, velocity, temperature, entropy and volumetric-heat-release rate. The 
bulk gas Is characterized by two equations: 
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where 
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Equations (11 to (5) are linearized by assuming each of the variables $, 
u, x, s, and q consists of a steady component that Is a function of axial 
position designated by the subscript ( )q and a small perturbation that Is 
a function of axial position and time designated by ( )i. The steady com- 
ponents are assumed to change gradually. Consequently, In small regions they 
can be considered to be constant. The steady state components satisfy the 
following zeroth-order system of equations: 
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The solution to equations (8) to (10) was obtained by direct Integration with 
the Runge-Kutta method (ref. 40) after rewriting the system of equations as 
follows: 
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In performing the Integration the steady state Inlet boundary conditions are 
selected to match the measured Inlet mass flow, pressure and temperature con- 
ditions and a volumetric-heat-source distribution Is selected so that the cal- 
culated temperatures In the long duct are similar to the measured temperatures. 


The perturbed components satisfy the following first-order system of 
equations: 
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^In order to solve equations (16) to (18) the periodic solutions are 
assumed to have the form; 


p^e.x) = P-| ( x) e~ 1<oe 

(19) 
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Note that the fluctuating volumetrlc-heat-f lux rate Is assumed to be a constant 
Independent of frequency and to be located on a disk at x s . Substituting 
equations (19) to (22) Into equations (16) to (18) yields the following equa- 
tion In matrix form for the spatial distribution of the perturbed quantities: 
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where T Is the transpose. 


Nonsource Region 


In the region where the unsteady heat release rate' Is zero, the resulting 
homogeneous differential equation Is 
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For the general case, an approximate solution to equation (30) can be 
obtained In terms of a transmission matrix by dividing the nonsource region 
Into a number of regions or subsections where steady-state quantities such as 
the Mach number and heat-release rate In each region can be assumed constant 
and the area variation In each region can be approximated by an exponential 
area variation. Note that In applying these equations to the particular test 
configuration used the area Is constant and the area variation Is zero. 

With these assumptions the B matrix Is Independent of x In the evalua- 
tion region. Thus the solution to equation (30) for the region Is 
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where ax Is the length of the region and exp(BjAxj) Is the matrix 
exponential. The matrix exp(BjAxj) Is the transmission matrix of the 
region ( T) j so that 
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Consequently, In a region without sources the transmission matrix Is found 
from 
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and N Is the number of regions. While the series given In equation (42) 
defines the matrix exponential. It was not used for computer calculation since 
the convergence of the series may be slow. The matrix exponential was evalu- 
ated from 
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and P Is determined to make the D matrix Is diagonal. The P matrix con- 
sists of the eigenvectors of sA. 


The P matrix used has the following elements for J = 1,2, and 3 
P(l.J) = (**^33 - x j) (^22 " X p “^2*^23 
P(2,J) = - |W 21 *^33 - Vj) -^ 31 ^ 23 j 

P(3,j) =-^21^32 "^31 ^22 ~ X j^ 

The values of x-j , \2» an< ^ x 3 are the eigenvalues of the matrix. 
These values are the roots of 
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Nozzle Boundary Condition 
The duct exit Impedance used Is 


pM(j) _ Z 

Vo u(1w) = Vo 


= 1 

EXIT 


V V 



(53) 


(54) 


where 

p * 2 k r (55) 

J-j ( p) Is the Bessel function of the first order and first kind, S(p) Is 
the Struve function of the first order and first kind and r Is the nozzle 
exit radius. This Is the Impedance of a circular duct with no flange having 
an open end fitted with an Infinite acoustically rigid flange (ref. 41). 


Source Region 

The source region solution Is obtained by Integrating equation (23) over 
a small Interval, c. 
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At the combustor Inlet the Impedance Is written In terms of a reflection 
coefficient: 
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Also at the Inlet the gas Is assumed to behave adlabatlcal ly 
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The resulting linear system of equations consists of three equations and has 
three unknowns. The unknown quantities are the pressure perturbations at the 
combustor Inlet and the duct exit and the pressure-temperature perturbation 
relationship at the duct exit. Consequently, for a given heat-release-rate 
perturbation the quantities can be determined as a function of frequency. 


EXPERIMENTAL INVESTIGATION 

The temperature and pressure measurements were made In the combustion rig 
shown In figure 1, which consisted of a J-47 can combustor (fig. 1(b)) mounted 
In a 0.203-m-dlameter by 0.77-m-long test section and an attached 6.44-m-long 
by 0.203-m-dlameter stainless-steel duct. The fuel used was hydrogen. 

Pressure and temperature measurements at stations 3 and 4 are discussed. All 
tests were conducted at the outdoor acoustic arena shown In figure 2. Tests 
were run for a range of conditions. Presented are measurements made at two of 
the test conditions shown In table I. Note that for the higher air mass flow 
rate cases the long duct exit temperature Is greater than the Inlet temperature 
Indicating some burning Is occurring In the long duct. 
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Pressure Measurements 


The Internal pressure transducers used were conventional 0.635-cm micro- 
phones with pressure response cartridges. To avoid direct exposure of the 
microphone to the combustion gases, they were mounted outside the ducted com- 
bustion rig, and the fluctuating pressure In the rig was transmitted to the 
transducers by means of a "semi-infinite" acoustic waveguide. Details of these 
probes are given In references 42 and 43. 


Temperature Measurements 

The dynamic temperature fluctuations were measured with Chromel-Alumel 
thermocouples. The dual-thermocouple probe design shown In figure 3 was 
selected so that the response time of the thermocouple could be determined 
from the actual test data rather than from a separate precallbratlon. The 
calibration procedure discussed In reference 14 was used. The thermocouple 
wire diameters were 76.2 and 25.4 ym. The thermocouples were Inserted Into 
the duct by actuators after combustion started In order to prevent their 
destruction by the combustion startup transient. 

Mean temperatures were measured at the combustor Inlet and at the measure- 
ment stations In the long duct. No radial or axial temperature distribution 
measurements were made. 


Data Acquisition and Processing 

Signals from the pressure and temperature probes were simultaneously re- 
corded on magnetic tape for later processing. Signals were recorded for 2 
minutes at each operating condition. The microphones were calibrated each day 
with a standard plstonphone before data were recorded. The thermocouples were 
used simultaneously to measure mean and fluctuating temperature. The signal 
cross and auto spectra were obtained- by off-line processing of the tape- 
recorded data on a two-channel fast-Fourler-transform digital signal processor 
with built-in analog-to-dlgltal converters and 120-dB/octave anti-aliasing 
filters. The resulting spectra were then transmitted from the processor to a 
central computer where an analyls of the spectra produced the thermocouple 
calibration. Then the calibration correction was Incorporated In a central 
computer program used to plot the data. 


COMPARISON OF THEORY AND MEASUREMENTS 

Measurements and theoretical calculations were made for the four test 
conditions shown In table I. However, similar test results were obtained at 
the low air mass flow rate (0.72 kg/s for Inlet temperatures of 822 K (Test 
115) and 1016 K (Test 116). Moreover, at the high air mass flow rate (1.54 
kg/s) similar test results were obtained for Inlet temperatures of 807 K (Test 
117) and 667 K (Test 118). Nevertheless, for each flow, the test results were 
different. Consequently, only one test condition Is shown at each flow rate. 
Since test l l 5 and 117 were done using similar duct Inlet temperatures, these 
results are presented and the results for tests 116 and 118 which were similar 
to those for tests 115 and 117, respectively, are not shown. For each case 
the fuel was hydrogen gas. 
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At each test condition theoretical and experimental pressure auto and 
cross spectra are compared to check the acoustic characteristics of a particu- 
lar model. Next, at each test condition theoretical and experimental tempera- 
turesauto and cross spectra are compared to check the thermal characteristics 
of a particular model. Last, at each test condition theoretical and experi- 
mental pressure-temperature cross spectra are compared to check the source 
region specifications of a model. Consequently, the following six types of 
theoretical and experimental results are compared for each operating condition: 

(1) Pressure auto spectra at the duct Inlet; 

(2) Pressure auto spectra at the duct exit; 

(3) Pressure cross spectra between the duct Inlet' and the exit; 

(4) Temperature cross spectra between the duct Inlet and the exit; 

(5) Pressure-temperature cross spectra at the duct Inlet; 

(6) Pressure-temperature cross spectra at the duct exit. 

Furthermore, for each cross spectrum the corresponding coherence function 
Is also presented. Note that a coherence function smaller than unity Implies 
one or more of the following conditions exists: 

(1) Extraneous noise Is present In the measurements. 

(2) The system relating the Input to the output Is nonlinear. 

(3) The output Is due to more than one Input. 

As previously discussed the zeroth-order system of equations (eqs. (11) 
to (13)) Is solved by Runge-Kutta Integration In order to establish the mean 
pressure, temperature, and velocity In the duct. Then, the transfer matrix 
equations can be used to determine the fluctuating quantities. For all the 
theoretical calculations the following steady-state mean volumetrlc-heat-f lux- 
distribution source equation was assumed and Is used In equations (11) and 
( 12 ): 



e~ nX - 10 3 (t - 288) 

g 


(65) 


The first term represents the heat Input In the combustor. The second term 
represents heat addition downstream of the combustor. The third term repre- 
sents heat loss through the wall of the long duct. The axial temperature dis- 
tribution Is not known; consequently, the parameters Qg, Q-j , <t, and n In 
equation (65) are arbitrary. They are established by an Iteration procedure 
which for each flow condition and each estimate of the parameters compares the 
calculated and measured temperature at the Inlet and exit of the long duct. 
Calculations were made for cases of gradual, sharp, and small temperature rise 
In the combustor to the duct Inlet temperature. For most comparisons of theo- 
retical and experimental results the theoretical calculations for only the 
gradual and sharp temperature rise models are presented. For each temperature 
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rise model and at each test condition the values of Qg and Q-| In 
equation (65) were adjusted until the calculated Inlet and exit temperatures 
were similar to those measured using a combustor Inlet temperature, tginj, of 
288 K, and a values of 0.5 and 0.1. The third theoretical calculation used a 
combustor Inlet temperature much closer to the duct Inlet temperature with a 
a value of 0.5. This case represents a situation where mixing of combustion 
gases and Inlet air yields a more uniform temperature In the combustor. Again 
Qg and Qi were adjusted so that the measured and calculated duct Inlet 
and exit temperatures were similar. For all cases x 0 = 0.6 and n ■ 0.320 
and the measured mass flow rate was used. The parameters and resulting tem- 
peratures for each theoretical calculation are given In table II. 

These three cases were selected to create examples having a sharp, a 
gradual, and a small temperature change In the combustor region while keeping 
the same temperature distribution In the long duct. Figure 4 shows the three 
resulting axial temperature profiles for test case 115. With tgi^ = 288 K and 
<s = 0.5 the temperature gradually Increases to t^LET- For tgjig = 288 K and 
a a 0.1 the temperature Increases abruptly to Tjnlet- For tg^ - INLET 
and a = 0.5, temperature changes In the combustor section are small. 

The combustor Inlet reflection factor for each case was assumed to be a 
constant given by 

R CIN = (0 ' 78 + 1 °- 550 > (66) 

This value was selected after many trial calculations. The heat perturbation 
source was assumed to be near the combustor Inlet. 

All calculations used the same value for the fluctatlng heat-release rate 
since the major objective was to compare’ measured and calculated pressure- 
temperature cross spectra. Consequently, the pressure auto-spectra are shown 
without levels. The level of the spectrum changes linearly with the value of 
the square of the f luctuatlng-heat-release rate (eq. (26)). However, the 
selecting of the proper f luctuatlng-heat-release rate for a significant number 
of cases, finding the relation between the fluctuatlng-and steady-state-heat- 
release rates, and studying the variation of the f luctuatlng-heat-release rate 
with frequency were tasks beyond the scope of the research program. 

First, acoustic and thermal theoretical and experimental results are pre- 
sented to demonstrate that the the model can be used to calculate these quanti- 
ties and to show that the source region model only has a secondary effect on 
these quantities. Then the pressure-temperature cross spectra results are 
presented. These show clearly the Importance of the selection of of the source 
region model. 


Pressure Auto Spectra 

Measured and calculated pressure auto spectra at the duct Inlet and exit 
for low velocity test case 115 and high velocity test case 117 are shown in 
figures 5 and 6. All the data presented were measured using a constant- 
bandwidth of 0.8 Hz. The plots show cancellations and reinforcements due to a 
longitudinal standing wave pattern. 
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Comparing calculated Inlet and exit pressure spectra for tciN = 288 K 
and a = 0.5 (figs. 5(b) and 6(b)) with those for tdN * 288 K and a = 0.1 
(figs. 5(c) and 6(c)) shows that the abrupt change In temperature associated 
wlthx a = 0.1 produces a spectrum that Is more Irregular than the spectrum 
associated with a = 0.5. Using a uniform temperature profile (figs. 5(d) and 
6(d)) yields a spectrum that Is more damped than the spectrum associated with 
either a = 0.5 or a = 0.1. Consequently, the small temperature gradient 
assumption does not appear appropriate. 

. For all cases the measured and calculated duct Inlet pressure auto spectra 
have cancellations at Identical frequencies and the shapes of the peak at the 
first reinforcement are similar. Also, between 40 Hz and 200 Hz the shape of 
the measured and calculated peaks at the reinforcements are not similar. How- 
ever, above 200 Hz the measured and calculated Inlet pressure spectra again are 
similar. 

The measured pressure spectra at the duct exit for low velocity test cases 
115 (fig. 5(a)(11)) and 116 (not shown) have similar combustion excited peaks 
near 120Hz which are missing In high velocity test cases 117 (fig. 6(a)(11)) 
and 118 (not shown). These peaks are thought to be due a self-excited combus- 
tion noise mechanism. Except for the small combustion temperature rise case, 
the calculated pressure spectra at the duct exit for low velocity test cases 
115 (fig. 5) and 116 (not shown) have their largest peaks at frequencies below 
120 Hz while for the high velocity cases the largest peaks are at higher 
frequencies. 

Examining the measured exit pressure spectrum shows that peaks are damped 
below 80 Hz for low velocity test case 115 (fig. 5(a)(11)) and below 160 Hz 
for high velocity test case 117 (fig. 6(a)(11)). The calculated exit pressure 
spectra do not exhibit this behavior. Finding the physical mechanism respon- 
sible for this behavior and removing this discrepancy between the measured and 
calculated exit pressure spectra Is one of the most Important tasks remaining 
In modeling this combustion noise data. 

The calculated and measured peaks and dips In the exit pressure spectra 
do not match for the low velocity cases. However, they are In agreement for 
the high velocity cases. This type of discrepancy may be due to use of an 
Incorrect temperature profile In the long duct. Only measurements at the Inlet 
and exit of the long duct were used to determine the temperature profile and 
this may not be sufficient when combustion occurs In the long duct. 

While these results show room for Improvement In calculating the duct 
exit pressure spectra, the model does a fair Job at calculating the duct Inlet 
pressure spectra. This aspect of the calculations Is Important In the present 
paper since the objective Is to compare measured and calculated pressure- 
temperature cross spectra at the duct Inlet and to find a single model that 
can be used to calculate both the pressure and temperature dynamics of the 
system. 


Pressure Cross Spectra 


Measured and calculated pressure cross spectra between the duct Inlet and 
exit are shown In figures 7 and 8 for low and high velocities (Test cases 115 
and 117), respectively. The measured and calculated phase angle plots show 
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phase angle transitions at Identical frequencies. Also, the slopes of the 
measured and calculated phase angle plots are similar for each case. In addi- 
tion, the shapes of the first peak In the measured and calculated cross spectra 
are similar. 

Again, the shapes of the measured and calculated peaks at the reinforce- 
ments are not similar In the frequency range from 40 to 200 Hz and are similar 
above 200 Hz. Comparing calculations for tciN = 288 K and a = 0.5 with 
those for tcmj = 288 K and a = 0.1 shows again that the a » 0.1 case 
produces a more Irregular spectrum than the spectrum calculated using a = 0.5. 
Although It Is not shown, as with the auto spectra, the uniform (small temper- 
ature rise) combustion temperature profile assumption yields a more damped 
spectrum than any other case. 


Temperature-temperature Coherence and Cross Spectra 

The coherence between temperatures measured near the duct Inlet and exit 
Is shown In figures 9 and 10 for low velocity test case 115 and high velocity 
test case 117, respectively. The coherence Is near 0.1 at the low frequencies. 
Note that spikes occur at 60 Hz Indicating some electrical power line noise Is 
present In the thermocouple signals. 

Measured and calculated temperature-temperature cross spectra between the 
duct Inlet and the exit are shown In figures 11 and 12 for test case 115 and 
117, respectively. The reinforcements and cancellations predicted In the tem- 
perature cross spectra magnitude are not observable In the measurements. How- 
ever, at frequencies up to 40 Hz for the low air mass flow rate case and 80 Hz 
for the high air mass flow rate the phase angles of both the measured and cal- 
culated results change linearly with frequency with a slope proportional to 
the time delay due to convection of temperature disturbances In the duct. 

At frequencies above 120 Hz the magnitude Is flat and the phase angle Is 
random because the response of the thermocouples Is reduced so much that the 
signals are at the level of the system electronic noise. In the corresponding 
uncompensated cross spectrum (not shown) the magnitude at 120 Hz Is 40 dB down 
from Its peak value. 

The lack of signal causes the coherence shown In figures 10 and 11 to be 
very small above 120 Hz. The measured phase angles shown In figures 11(a) and 
12(a) are more erratic when the coherence Is small. Note, however, that at 
coherence values greater than 0.01 good phase angle values are obtained. 


Pressure Temperature Coherence and Cross Spectra 

The pressure-temperature coherence at the duct Inlet Is shown In 
figures 13 and 14 for test cases 115 and 117, respectively. As just noted, at 
low frequencies the coherence Is large enough (above 0.01) so that the cross 
spectra phase angle measurements will be valid. This can also be Inferred 
from a comparison of the phase angle plots shown In figures 15 and 16 and the 
corresponding coherence plots. Again, the pressure-temperature cross spectra 
phase angle values are more erratic when the coherence Is small. 
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Measured and calculated pressure temperature cross spectra at the duct 
Inlet are shown In figures 15 and 16 for test case 115 and 117, respectively. 
The locations of the cancellations In the measured and calculated cross spectra 
magnitude plots are similar. However, the shapes of the peaks near the rein- 
forcements are not similar. The cross spectra phase angle plots are similar 
below 160 Hz. for the low air mass flow rate case and below 200 Hz for the 
high air mass flow rate case. 

Pressure temperature coherence plots at the duct exit are shown In 
figures 17 and 18 for test case 115 and 117, respectively. At most points the 
coherence Is so low that the cross spectra phase angle measurements are not 
reliable. 

Measured and calculated pressure temperature cross- spectra at the duct 
exit are shown In figures 19 and 20 for test cases 115 and 117, respectively. 

At most points the correspondence between the calculated and measured cross 
spectra shown In figures 19 and 20 Is negligible. 

The phase angle relationship Is controlled by the convective velocity and 
the speed of sound along the duct. Both quantities are temperature dependent 
and thus the convective velocity and speed of sound are functions of the tem- 
perature distribution. The theoretical pressure-temperature cross spectra 
near the duct Inlet obtained for a gradual temperature change rather than an 
abrupt temperature change show the best agreement with the measurements at the 
higher air mass flow rates (fig. 16(b)). However, at the lower air mass flow 
rates the uniform (small temperature rise) combustor temperature profile model 
calculations yields the best agreement (fig. 15(d)). Since these results In- 
dicate that the small temperature rise model Is best for the pressure tempera- 
ture cross spectra at the low flow rate, while the pressure auto-spectra 
calculations (fig. 6(d)) Indicate It Is not, the overall physical model Is 
least successful for the low velocity cases. 

The theoretical calculations show that at a given air mass flow rate the 
amount and slope of the temperature change In the source region determines the 
pressure-temperature phase angle variation and Is more Important than the Inlet 
and exit temperatures In determining the phase angle curve shape. This 
behavior Is In good agreement with the measured pressure-temperature phase 
angle plots. 

These results show that the model, which explicitly links dynamic- 
volumetric-heat-source fluctuations to pressure and temperature fluctuations. 

Is In good agreement with the measurements made at the hlger flow rate. The 
correlation of combustion noise measurements with a single parameter related 
to the heat-release rate that can be calculated from the steady state operating 
conditions and the geometry, thus In effect, links the steady-state volumetric 
heat-release rate and the dynamic-heat-source fluctuations. The pressure- 
temperature cross spectra shown In figures 15 and 16 show, on a theoretical 
basis rather than from a statistical one, that a definite physical relationship 
between the pressure and the heat source exists. These results also Indicate 
that the physical mechanism of Interest Involves the steady-state heat flux 
and the dynamic fluctuating heat Input. In other words, the many correlations 
that have been developed that link combustion noise and heat flux are actually 
relations between the steady-state volumetric heat-flux and the dynamic volu- 
metric f luctuatlng-heat Input. 
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CONCLUSIONS 


A one dimensional treatment of unsteady pressure and temperature dlsturb- 
ancesgeneratlon and propagation In a combustion duct has been described. 
Equations that Include a dynamic-volumetric-disk-heat-flux-rate source and a 
spatially distributed steady-state volumetric-heat-flux rate are derived. 

The equations were solved by matrix methods. Good agreement was obtained 
between the measurements and the calculations for the higher of the two air 
mass flow rates used (1.57 kg/s). The results show that combustion noise Is 
caused by fluctuating volumetric-heat-flux rate. Improved results might be 
obtained If more were known about the spatial distribution of the steady-state 
volumetric-heat-release rate, the acoustic Impedance boundary conditions at 
the combustor Inlet and duct exit, self-excited combustion noise, and the 
physical mechanism that damps the duct exit pressure spectrum at the low 
frequencies. 

The use of the dynamic themocouple data does provide enough new Informa- 
tion to suggest that the model used gives predictions that are accurate enough 
for practical purposes, provided the frequency of Interest Is low enough so 
that only the plane wave propagates. 

The results show that correlations between parameters related to the 
steady state volumetric-heat-release rate and combustion noise are probably 
more correctly described as correlations between the steady-state volumetric- 
heat-flux rate and the dynamic f luctuatlng-volumetrlc-heat-f lux rate. The 
experimental and theoretical results shown here demonstrate the physical link 
between the fluctuating volumetric-heat-release rate and combustion acoustics. 
Consequently, the combustion noise mechanism that should be Investigated to 
better understand combustion noise Involves the steady and dynamic volumetric- 
heat-flux rate In various combustors. 
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TABLE I. - TEST CONDITIONS 


Test condition 


Test 



115 

116 

117 

118 

Air mass flow rate, kg/s 

0.72 


1 .54 

1.55 

Hydrogen mass flow rate, kg/s 

0.0031 




Fuel-air ratio 
Duct temperature, K: 

0.0043 


0.0045 


f Inlet 

822 

1016 

807 

667 

T ex1t 

792 

954 


765 


TABLE II. - THEORETICAL CALCULATIONS 


Condition 






Test 

No. 








115 



116 



117 



118 






Hass flow rate, W 

. kg/s 






0.72 



0.73 



1.54 



1 .55 


a 

0.5 

0.5 

0.1 

0.5 

0.5 

0.1 

0.5 

0.5 

0.1 

0.5 

0.5 

0.1 

Combustor Inlet temperature, 
T c,1n> K 

Volumetric-heat- release-rate 
coefficients : 

288 

795 

288 

288 

995 

288 

288 

780 

288 

288 

620 

288 

Q 0 x 106 

12.8 

0.9 

12.8 

18.3 

0.9 

18.4 

24 

0.1 

24.6 

16.7 

0.1 

17.1 

QtxIO 6 

Duct temperature, K: 

1 . 3 

1 .3 

-1.35- 

1.3 

1.3 

1.4 j 

1 

5.2 

5.0- 

5.4 

4.8 

5.0 

4.95 

Tlnlet 

826 

826 

822 

1018 

1013 

1016 

800 

808 

808 

661 

667 

667 

T ex1t 

795 

793 

792 

954 

946 

954 

892 

904 

902 

758 

766 

763 
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(a) RIG SCHEMATIC. 



(b) INJECTION SCHEME. 


Figure 1 . - Schematic of ducted combustion system rig. 




Figure 2. - NASA Lewis combustion acoustics facility. 
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Figure 3 0 - Dual-thermocouple probe 
design. 
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(ii) Pressure auto spectrum P 4 P 4 ' near 
duct exit. 

la) Measured pressure auto spectra for 
bandwidth of 0 . 8 Hz. 

Figure 5. - Pressure auto spectrum for 
test 115. 
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(i) Pressure auto spectrum P 3 P 3 * near 
duct inlet. 



FREQUENCY, F, Hz 


(ii) Pressure auto spectrum p 4 p 4 * near 
duct exit. 

(b) Calculated auto spectra - 288 K, 

0 - 0. 5). 


Figure 5. - Continued. 
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(i) Pressure auto spectrum P 3 P 3 * near 
duct inlet. 



FREQUENCY, F, Hz 


(ii) Pressure auto spectrum p^p/ near 
duct exit. 

(c) Calculated spectra (tQ| N 0 288 K, 

o=0.1.) 


Figure 5. - Continued. 


PRESSURE LEVEL, dB 



(i) Pressure auto spectrum P 3 P 3 * near 
duct inlet. 



(ii) Pressure auto spectrum p 4 p 4 * near 
duct exit 


(d) Calculated spectrum (t£| N = 795 K, 
o = 0.5). 


Figured - Concluded. 


PRESSURE LEVEL, dB (re. 20 jjPa) 



(i) Pressure auto spectrum P 3 P 3 * near 
duct inlet. 



(ii) Pressure auto spectrum p^p/ near 
duct exit. 

(a) Measured pressure auto spectra for 
bandwidth of 0. 08 Hz. 

Figure 6 . - Pressure auto spectra for 
test 117. 
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(i) Pressure auto spectrum P 3 P 3 * near 
duct inlet. 



(ii) Pressure auto spectrum p 4 p 4 * near 
duct exit. 

lb) Calculated auto spectra (tQ| N = 288 K, 

0 = 0. 5). 


Figure 6 . - Continued. 
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(i) Pressure auto spectrum P 3 P 3 * near 
duct inlet. 
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(ii) Pressure auto spectrum p 4 p 4 * 
near duct exit. 

(c) Calculated auto spectra (tQ^ - 288 K, 

0 - 0.1). 


Figure 6 . - Continued. 
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(i) Pressure auto spectrum P 3 P 3 * near 
duct inlet 



FREQUENCY, F, Hz 


(ii) Pressure auto spectrum P 4 P 4 * near 
duct exit. ...... 

(d) Calculated auto spectra (tQ JN = 780 K, 
o = 0.5). 

Figure 6 . - Concluded. 
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(a) Measured pressure cross spectrum 
P 4 P 3 * for bandwidth 0 . 8 Hz. 

Figure 7. - Pressure cross spectra between 
station 3 and 4 for test 115. 
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(b) Calculated pressure cross spectra 
P4P3* {t CIN = 288 K . 0 =0.5) 


Figure 7. - Continued. 



(c) Calculated pressure cross spectra 
P 4 P 3 * ^j|\j = 288 K, 0 = 0. 1). 

Figure 7. - Concluded. 
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(a) Measured pressure cross spectrum 
P4P3* bandwidth of 0. 8 Hz. 

Figure 8. - Pressure cross' spectra between 
station 3 and 4 for test 117. 
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lb) Calculated pressure cross spectra 
P4P3* ^CIN = 0 = 0. 5). 

Figure 8. - Continued. 
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(c) Calculated pressure cross spectra 
P4P3* (tcii\| = 288 K, 0 = 0. 5). 


Figure 8. - Concluded. 



FREQUENCY, Hz 

Figure 9. - Coherence between duct in- 
let and exit temperatures 1 1 4 I 2 / 

1 1 4 1 2 1 to 1 2 , test no. 115 (bandwidth = 
0.8 Hz). 


MAGNITUDE, dB 

PHASE, deg (re. 4.4|iK) COHERENCE 



FREQUENCY, F, Hz 


Figure 10. - Coherence between duct in- 
let and exit temperature I I 2 / 1 1 4 1 2 
1 1} I 2 for test 117 (bandwidth = 0. 5 h3. 
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(a) Measured temperature cross spectra 
with a bandwidth of 0. 8 Hz. 


Figure 11. - Temperature cross spectra 
t^* between duct inlet and duct exit 
for test 115. 
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(b) Calculated temperature cross spectrum 
tt c|N = 288K, 0 = 0.5). 


Figure 11. - Continued. 
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(a) Measured temperature cross spectra 
with bandwidth of 0. 8 Hz. 

Figure 12. - Temperature cross spectra 
t^* between duct inlet and duct exit 
for test 117. 
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(b) Calculated temperature cross spectrum 
ttciN = K, o = 0. 5) 


Figure 12. - Continued. 
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Figure 13. - Coherence between pressure 
and temperature at duct inlet l p 3 t 3 * I 2 
/ 1 p 3 1 2 1 1 3 1 2 for test 115(bandwidth = 

0. 8 Hz). 
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Figure 14. - Coherence between pressure 
and temperature at duct inlet I 1 2 

/ I P 3 I 2 1 13 I 2 for test 117 (bandwidth ■ 
0. 8 Hz). 
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(a) Measured pressure temperature cross 
spectrum with a bandwidth of 0. 8 Hz. 

Figure 15. - Pressure temperature cross 
spectra near duct inlet poto# for test 
115. 
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(c) Calculated pressure temperature 
cross spectrum (t^| N = 288 K, o - 0. 1). 

Figure 15. - Continued. 



_U -—J I ' - • i ■ ' A 1 • I • ' l_ ' . .1 


0 40 80 120 160 200 240 280 320 



PHASE ANGLE, deg MAGNITUDE, dB 



0 40 80 120 160 200 240 280 320 


FREQUENCY, F, Hz 

(d) Calculated pressure temperature 
cross spectrum (t^ = 795 K, o = 0. 5). 

Figure 15. - Concluded. 
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(a) Measured pressure temperature 
cross spectrum with a bandwidth of 
0.8 Hz. 

Figure 16. - Pressure temperature cross 
spectra near duct inlet P 3 t 3 * for test 117. 
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(b) Calculated pressure temperature 
cross spectrum (t^ = 288 K, o = 0. 5). 
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(c) Calculated pressure temperature 
cross spectrum (t^ = 288 K, o = 0. 1). 

Figure 16. - Concluded. 
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Figure 17. - Coherence between pressure 
and temperature near duct exit, 1 p 4 t/ 1 2 
/ i P 4 i 2 1 1 4 1 2 for test 115 (bandwidth = 

0.8 Hz). 
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Figure 18. - Coherence between pressure 
and temperature near duct exit, I I 2 
/ 1 p 4 1 2 1 1 4 1 2 for test 1 17 (bandwidth = 
0.8 Hz). 
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(a) Measured pressure temperature cross 
spectrum with a bandwidth of 0. 8 Hz. 

Figure 19. - Pressure temperature cross 
spectra near duct exit p^* for test 115. 
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(b) Calculated pressure temperature cross 
spectrum (t^ = 288 K, o = 0. 5). 

Figure 19. - Continued. 
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(c) Calculated pressure temperature cross 
spectrum (Tqj N = 288 K, o = 0. 1). 

Figure 19. - Concluded. 
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(a) Measured pressure temperature cross 
spectrum p^* for a bandwidth of 0.8 Hz. 

Figure 20. - Pressure temperature cross 
spectra near duct exit p^* for test 117. 


( 



0 40 80 120 160 200 240 280 320 


FREQUENCY, F, Hz 

(b) Calculated pressure temperature 
cross spectrum (tc| N = 288, o = 0. 5). 

Figure 20. - Continued. 
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(c) Calculated pressure temperature cross 
spectrum {tQj N = 288 K, a - 0. 1). 

Figure 20. - Concluded. 
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